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The validity of the retracing approximation in the semiclassical quantization of Andreev billiards 
is investigated. The exact energy spectrum and the eigenstates of normal-conducting, ballistic 
quantum dots in contact with a superconductor are calculated by solving the Bogoliubov-de Gennes 
equation and compared with the semiclassical Bohr-Sommerfeld quantization for periodic orbits 
which result from Andreev reflections. We find deviations that are due to the assumption of exact 
retracing electron-hole orbits rather than the semiclassical approximation, as a concurrently per- 
formed Einstein-Brillouin-Keller quantization demonstrates. We identify three different mechanisms 
producing non-retracing orbits which are directly identified through differences between electron and 
hole wave functions. 
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I. INTRODUCTION 

The dynamics of hybrid systems consisting of a normal- 
conducting (N) quantum dot and a superconducting (S) 
waveguide have recently raised much experimental^ 2 ^ 
and theoretical 4 i 5 i 6 i 7 i 8 i 9 i 1Q i n i 12 interest. This is due to 
the unusual and sometimes counterintuitive properties 
of the interface of normal- and superconducting regions. 
An electron moving inside the N region will be co- 
herently scattered into a hole upon contact with the 
superconductor [see Fig. [Ha)]™ This phenomenon is 
commonly referred to as Andreev reflection™ In the 
frequently used and remarkably successful semiclassical 
Bohr-Sommerfeld (BS) methocU 5 ^^ it is assumed that 
the Andreev reflection is perfect, i.e. that the path of the 
backscattered hole will exactly trace that of the incoming 
electron (v e = —Vh) [see Fig. QJb)]. The consequences 
of this assumption are profound: all trajectories ema- 
nating from the SN-interface are strictly periodic. The 
classical dynamics of the combined SN-system become 
entirely regular, even and in particular when the nor- 
mal conducting cavity would feature hard chaos™ Un- 
like chaotic or integrable systems, periodic orbits are no 
longer isolated but form continuous manifolds that dom- 
inate the classical phase space and, in turn, the den- 
sity of states (DOS). 19 ! 20 ! 21 The BS approach 15 ! 16 ! 22 to 
the DOS of an Andreev billiard relies on three assump- 
tions: exact retracing of electron-hole trajectories (re- 
ferred to as assumption Al in the following), the absence 
of any quasi-periodic orbits other than the ones caused 
by Andreev reflection (assumption A2), and the applica- 
bility of semiclassical approximations (assumption A3), 
i.e. Xd <C \/A, where is the de Broglie wavelength 
and A is the area of the N billiard. 

The BS approximation has been found to be supris- 
ingly accurate^SS even in the presence of a magnetic 
fielo™ or a soft wall potential™ However, some difficul- 
ties have emerged: for example, the excitation gap in 
the DOS for billiards with a chaotic TV cavity cannot be 
properly accounted f or ^iSSiS& A direct inspection of the 




FIG. 1: (a) Andreev billiard with a rectangular normal- 
conducting (N) region. The superconducting (S) lead is as- 
sumed to be half-infinite (shaded area). The dashed and dot- 
ted lines depict an "almost retracing" electron-hole orbit cre- 
ated by Andreev reflection at the SN-interface. For better 
visibility the starting point of the orbit is marked by a dot. 
(b) perfect Andreev reflection, i.e. e = Oh- (c) imperfect 



Andreev reflection, e 7^ Oh, ky 



electron- and hole-components of the wave functions 23,27 
revealed markedly different patterns, clearly signaling a 
breakdown of the retracing approximation. 

In the present communication, we investigate the 
role of non-retracing orbits in Andreev billiards. Non- 
retracing electron-hole orbits leave their mark on the An- 
dreev wave functions which we determine by solving the 
Bogoliubov-de Gennes equation with the help of the mod- 
ular recursive Green's function method (MRGM)™iSi 
We find that Andreev states which correspond to non- 
retracing orbits break the close correspondence between 
the electron and hole wave function patterns. The anal- 
ysis of the eigenstates allows us to check the merits and 
limitations of the BS approximation against the exact 
quantum mechanical approach. We identify different 
sources for the breakdown and explore separately each 
of the assumptions Al to A3 listed above. 

This paper is organized as follows. In Sec. |H] the 
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method for solving the quantum mechanical eigenvalue 
problem as well as the Bohr-Sommerfeld approach to An- 
dreev billiards is briefly reviewed. In section ITTTI we iden- 
tify different mechanisms leading to non-retracing orbits 
in Andreev billiards. The paper concludes with a short 
summary in Sec. IIVI 



does not rely on the usual Andreev appoximation, namely 
A <C Ep, and the assumption of approximately perpen- 
dicular angles of incidence on the SN- inter face Jiii This 
feature will be important when we consider the effect of 
imperfect Andreev reflections on the quantum and semi- 
classical spectra. 



II. MODEL SYSTEM 

We consider as a model system a rectangular nor- 
mal conducting cavity of width D, length L, and area 
A = LD which is attached to a half-infinite lead of width 
W [see Fig. Ufa)]. For technical reasons we focus on an 
A^-cavity that features classically regular motion. Since 
within the framework of perfect Andreev reflections com- 
bined SN- Andreev systems become regular and, in fact, 
periodic irrespective of the underlying regular or chaotic 
dynamics of the iV-cavity, we expect many of our results 
to be valid for arbitrary A"-cavities. As we consider only 
one superconducting lead, we choose the gap parameter 
A to be real. The superconducting coherence length (£) 
is assumed to be small compared to the linear dimension 
of the rectangular cavity (~\/]4). Under that assumption 
we may use a step function model 7 at the SN- interface, 
A = |Ao| 0(xns —x). This SN-interface, which is located 
at xns — 0, is assumed to be ideal, meaning that there is 
no tunnel barrier and no mismatch between Fermi ener- 
gies and effective masses in the N- and S-region. An exact 
treatment of the Bogoliubov-de Gennes equation reveals 
that this results in a probability for Andreev reflection 
close to unity (up to corrections of the order A 2 / E^) in 
the energy range we consider, < e < A = 0.021^^*^ 
Furthermore, the phase coherence length of the scatter- 
ing process is assumed to be infinitely long, resulting in 
fully ballistic dynamics. 



A. Quantum mechanical solution 

We describe the S-N hybrid system by the Bogoliubov- 
de Gennes (BdG) equation 30 
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Ho = p 2 /(2m e ff) + V(r) — Ep is the single particle 
Hamiltonian. The electron (hole) quasiparticle wave- 
functions are denoted by u (i?), e is the excitation en- 
ergy of the electron (hole) above (below) the Fermi en- 
ergy Ep. Throughout this publication we consider only 
bound Andreev states for which the solutions of in 
the superconducting region are exponentially decaying. 
To identify the bound Andreev states, we use a wave 
function matching technique^ in combination with the 
Modular Recursive Green's function methodS2i2£ as out- 
lined in Ref. E3 

Note that the above procedure is exact in the sense that it 



B. Semiclassical analysis 

The intuitive picture of ideal Andreev reflections 
[Fig. Dp))] lends itself to a semiclassical description and, 
more specifically, to a periodic-orbit quantization. An 
electron approaches the SN-interface with wave number 

k e = ^/2(E F + s) (2) 

and angle e relative to the interface normal. The hole 
leaves the interface with wavenumber 

k h = ^2(E F - e) . (3) 

and a corresponding angle Oh- Since the component along 
the interface is exactly conserved because of translational 
invariance, fc 6j2/ = kh, y (or v eyV = —Vh, y ) the compo- 
nents normal to the interface will be, in general, different 
k e ,x 7^ kh, x -> leading to imperfect retracing [Fig. [He)]. 
Only in the limit e —> (or e/Ep —> 0) perfect Andreev 
reflection k eyX = kh, x ( or v e,x = —Vh,x) ensues. The 
retracing approximation consists now of the assumption 
Al of perfect reflection, v e , x = —Vh, x , for all e in the in- 
terval < e < A. The validity and breakdown of this 
assumption will be explored in the following. 

Under assumption Al all trajectories emanating from 
the SN-interface are strictly periodic. For A^-cavities fea- 
turing hard chaos, every trajectory will eventually hit 
the SN-interface, thus yielding a globally periodic sys- 
tem. For Af-cavities with mixed or regular dynamics, 
certain regions of the phase space may remain decoupled 
from the interface and thus may feature both continu- 
ous manifolds of strictly periodic orbits and islands with 
quasi-periodic motion. Neglecting the latter contribution 
to the DOS implies assumption A2. 

The BS quantization of the continuum of perfectly pe- 
riodic orbits visiting the SN-interface yields for the action 
S of periodic orbits created by pairs of electron-hole tra- 
jectories (with wavenumbers fc e , kh), 16 

S = l(k e - k h ) = 2irn + 2 arccos(e/A). (4) 

The variable / stands for the length of an arbitrary path 
connecting one point at the SN-interface with another. 
The quantum number n = 0,1,2,... accounts for the 
quantized difference in action between the electron and 
the hole part of the periodic orbit. The second term on 
the right hand side of Eq. arccos(e/A), is the Maslov 
index 31 for one reflection at the SN-interface. Including 
this energy-dependent phase shift at the point of Andreev 
reflection extends^ the accuracy of the BS description to 
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energy ranges < e < A . An even more accurate semi- 
classical description of the Andreev reflection process will 
be discussed in a subsequent publication^ 

The semiclassical BS approximation for the state 
counting function can then be written as^fi 
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where l n (e) is defined by 
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The geometry of the cavity enters via the classical path 
length distribution P(l). The quantity P(l) dl is the clas- 
sical probability that an electron entering the N region 
at the SN-interface returns to the interface after a path 
length in the interval [/, I + dl]. 

In the special case that the Andreev system becomes 
separable, in particular for the geometry W = D [see the 
inset in Fig. 0a), and Fig. 0(a)], an alternative semiclas- 
sical approximation, the Einstein-Brillouin-Keller (EBK) 
quantization, becomes applicable. This alternative is of 
particular interest as it only relies on assumption A3, but 
does not involve the additional assumptions Al and A2, 
thus allowing to disentangle the validity of the semiclas- 
sical approximation from that of the retracing approx- 
imation. For an EBK quantization of the cavity with 
W = D the separability yields two quantization condi- 
tions, the first one of which pertains to the motion in x 
direction [compare with Eq. (|3J)], 



k*> h dx = 2L(k e x - *£) = 2 
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The second quantization condition in y direction reads 



kydy = 2k y D = 2mir . 



(8) 



The contour C x (C y ) stands for a path along a horizontal 
(vertical) line in the normal conducting billiard. Intro- 
ducing the short hand notation fc™(±£) = k^ with 
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the quantization condition for x [Eq. (J7J)] reads 



Tin + arccos | 



(^) — L i^m ~ km] 



(9) 



(10) 



This is a transcendental equation in e(n, ra), which al- 
lows us to calculate individual energy levels in a cavity 
with W = D semiclassically. Note that in contrast to 
this EBK result, the BS approach yields only an approx- 
imation to the smoothed state counting function N(s). 
We therefore transform the EBK results for e(n, ra) into 
a state counting function 
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and compare both semiclassical approximations with the 
quantum results. EBK quantization of Andreev billiards 
can be viewed as the analogue to the adiabatic quantiza- 
tion of smooth soft- wall chaotic billiards, for completely 
integrable systems As the EBK quantization invokes 
quasi-periodic rather than periodic motion no assump- 
tion of retracing is involved. 



III. VALIDITY OF THE RETRACING 
APPROXIMATION 

We will now explore the validity and breakdown of the 
retracing approximation and the resulting BS quantiza- 
tion by considering rectangular cavities (Fig.EJ) with dif- 
ferent ratios D/W and D/L. With this parameter space 
at our disposal we can probe and disentangle the validity 
of assumptions Al to A3. 



A. Quadratic N cavity 

We consider first the quadratic cavity with D = L. 
Here two cases arise: for D/W = 1, the Andreev billiard 
becomes separable while for D/W > 1 the Andreev sys- 
tem is non-separable even though the N cavity itself is 
separable. In the latter case the SN-boundary condition 
breaks the separability. The location of the placement 
of the SN-interface for D/W > 1 provides an additional 
degree of freedom for our investigation to be exploited 
below. A comparison between the BS approximation 
and the exact quantum mechanical calculation for the 
counting function N(s) is shown in Fig. [21 Overall, for 
D/W = 1 (a) and D/W = 1.25 (b), the BS approxima- 
tion performs quite well. In particular, it reproduces and 
intuitively explains the cusps seen in N(e). They are due 
to the sharp peaks in the pathlength spectrum. The po- 
sition of the cusps is predicted very accurately by the BS 
approximation with an error well below the mean level 
spacing of 0.09A. 

Assuming the validity of the retracing approximation, 
the semiclassical quantization condition [Eq. (JSJ] allows 
to map every excitation energy e onto a path length l n 
of a periodic Andreev orbit. Consider, e.g., the energy at 
the cusp marked / in Fig.[2{b) that corresponds to a path 
length Iq = 2L, which is the length of the shortest classi- 
cal Andreev-reflected orbit of the system [see P(l) in the 
inset of Fig. H(b)]: The electron leaves perpendicular to 
the SN-interface, is reflected at the opposing wall and re- 
turns to the SN-interface. Quantum mechanical wave 
functions evaluated at the cusp energy indeed feature 
a pronounced enhancement along the orbit "bundle"— 
with length I = 2L [see in Fig. [3(b)]. Wavefunctions 
with neigboring energy values below the cusp feature ad- 
ditional nodes in transverse direction [see Fig. 0(a)] . Note 
that also eigenstates near the cusp marked // in Fig.|3(c, 
d) correspond very nicely to the path-bundles of the ex- 
pected length. Consistent with the good agreement of 
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FIG. 2: (Color online) The quantum mechanical state count- 
ing function iVQM(^) (solid red staircase) and its semiclassical 
BS estimate Nbs(s) (dashed green line) for two quadratic 
cavities with different lead widths as shown in the insets 
[k F = 20.5tt/W, W = L(a), 0.8L(b)]. The top left insets 
show the classical path length distribution in units of the cav- 
ity length L. Solid triangles mark pronounced cusps in N(e) 
and their classical origin. The quantum number n (Eq. 0| 
increases from to 1 at /. Open triangles mark the energy 
positions of states whose wavefunctions are displayed in Fig.|3| 
Shading marks the regions where the BS approximation be- 
gins to deviate from quantum results. 



the BS quantization with the quantum calculations, the 
electron and hole wavefunction densities agree very well 
with each other. To the extent that bundles of classical 
trajectories cause the density enhancement in the wave- 
function, and bundles of hole- and particle-orbits agree 
with each other, this is to be expected. Conversely, a hall- 
mark of the breakdown of retracing are distinct density 
ditributions in the particle (u) and hole (v) wave func- 
tions, as has been recently observed £2*21 Looking more 
closely, we find discrepancies between the exact quantum 
mechanical calculations for Nqm(s) and its semiclassical 
counterpart Nbs(z), which are indicated by the shaded 
areas in Fig. |2{a) and (b). For the two systems con- 
sidered in Fig. [3 we note that the mismatch between the 
quantum and the semiclassical results is located at rather 
well-defined values of the excitation energy e. Note that 
the mismatch tends to occur at values of £ just above a 
cusp. 

Since pathlength and energy at fixed quantum number 





FIG. 3: (Color online) (a) — (d) shows the electron and hole 
wavefunctions |i/(x)| 2 and |v(x)| 2 at values of e indicated by 
open triangles in Fig. e/A is (a) 0.565, (b) 0.566, (c) 

0.33, (d) 0.355. The corresponding classical orbits are indi- 
cated on the right. 
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FIG. 4: Geometries with tunable boundaries: (a) stretched N- 
cavities with L ^> D but with D equal to the width of the SN- 
interface W = D. (b) stretched N-cavities with D ^> L but 
constant length L (lower boundary of S and N aligned), (c) 
as (b) but SN-interface at arbitrary position on the entrance 
side. 



n are inversely proportional to each other [Eq. 0] the 
energies above the cusp which is associated to quantum 
number n [e.g. the one marked / in Fig.EJb)] correspond 
to the longest orbits associated with quantum number 
n + l. The deficiencies of the BS approximation are ev- 
idently caused by contributions from long orbits. This 
observation is in line with the well-known failure of the 
BS approximation to yield the energy gap for chaotic N 
billiards ^i^iL also caused by long orbits. We also see a 
discrepancy above the cusp marked II in Fig.J^b) which 
turns out to be due to diffr active scattering at the corner. 
We will discuss diffractive effects in more detail below 
when we explore different mechanisms for the failure of 
the BS approximation by changing the geometry. 



B. Stretched separable billiard 

We consider now the separable case W = D but stretch 
the billiard to D/L <C 1 [see Fig. 0(a)] without, how- 
ever, changing the Fermi energy Ep of the system. The 
distance I travelled between two encounters with the SN- 
interface is, for the shortest orbits, at least 2L. The effect 
of imperfect retracing, originating from the difference in 
angles between e and Oh, is a divergence between elec- 
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tron and hole trajectories (ex /), which is amplified by 
large L. We introduce as a measure for the imperfect 
retracing the lateral displacement Sy on the interface be- 
tween the hole and the particle trajectory after one loop 
(see Fig. EJ, 



Sy = I |sin(0 e ) - sin(6^)| 



= I 



£sin( 



+ 0(£ 2 ), 



(12) 



where = e ^\ £ =o- For the particular geometry consid- 
ered [Fig. 0(a)], Sy increases linearly with both the tra- 
jectory length I and the excitation energy e. Note that 
the ratio e/E-p is, in general, much smaller than one, 
i.e. e < A <C E F . 

This classical scale for the mismatch between the par- 
ticle and hole orbits should be related to the quantum 
scale, i.e. the linear dimension of the wave packet esti- 
mated by the de Broglie wavelength Xp. We thus in- 
troduce r = \5y\/\p as order parameter for the error of 
retracing. For r«l the wave packet cannot resolve the 
mismatch and the BS approximation should work well. 
Conversely, as r reaches the order of unity, quantization 
based on the existence of a continuum of periodic orbits 
should fail. 

To probe the breakdown of the retracing approxima- 
tion quantitatively, we compare the semiclassical with 
the quantum density of states (DOS), which are obtained 
from the state counting function N(e) as 



Pbs,qm(^) 
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(13) 



since the DOS is more sensitive to errors than the 
(smoothed) spectral staircase. Results for two cavities 
with different values of L are shown in Fig. HJa) and 
(b). This example illustrates that the degree of agree- 
ment between pqm and pes is indeed controlled by r. 
The following trends can be observed: (i) Due to the 
comparatively shorter length I of orbits in (a), the over- 
all agreement there is better in (a) than in (b). (ii) The 
agreement in (b) deteriorates for larger values of s/A, 
since the mismatch in retracing increases with e. 




FIG. 5: Imperfect Andreev reflection in the fundamental 
and extended zone scheme. The returning hole hits the SN- 
interface a distance Sy apart from the starting point of the 
particle. 

As a measure for the average mismatch between pqm 
and pes in an Andreev billiard of given geometry we use 



the quantity Sy, which is Sy averaged over all e and all 
angles, 



1 r /2 2LA 
Sy = — — / (IB cos 1 de Sy = — — 

^ A J-tt/2 Jo 7rE~p 



(14) 



since for this particular geometry I = 2L/ cos . The error 
in the DOS is quantified in terms of the root mean square 
(RMS) deviation Sp between pBS and pqm, 



Sp- 



^ / de |p Q M - Pbs\ 2 ■ 



(15) 



As expected, as long as r = Sy/Xp <C 1, the retracing ap- 
proximation is sufficiently accurate to reproduce the DOS 
on the scale of the mean level spacing [see Fig.EJc)]. Note 
that resolution of individual levels is beyond the scope 
of the BS approximation. As discussed below, an EBK 
quantization is needed for an accurate quantization of in- 
dividual levels. As f approaches unity the oscillations in 
the quantum DOS cannot be resolved any longer. The er- 
ror Sp(Sy) appears to saturate, in agreement with the ob- 
servation that in this regime the main peaks in Fig. HJb) 
[corresponding to the main cusps in N(e)] remain well de- 
scribed by the retracing approximation. Identifying the 
latter with the shortest Andreev-reflected orbits of length 
~ 2L allows a simple explanation as to why even a drastic 
elongation of the cavity leaves these cusps well described 
by the retracing approximation: orbits of length ~ 2L 
correspond to strictly horizontally injected trajectories 
(with = 0). They do not accumulate any lateral dis- 
placement Sy irrespective of the length of the orbit and e. 
While the average value Sy is increasing with growing cav- 
ity elongation (due to orbits_with 0^0), the horizontal 
orbits limit the increase Sp(Sy), leading to the observed 
saturation in Fig.EJc) for Sy > Af- 

Additional evidence that the discrepancy between the 
quantum DOS and the BS spectrum is due to the re- 
tracing assumption (Al) and not due to the failure of 
semiclassics itself (A3) can be gained from an EBK quan- 
tization. Note that the Andreev system [Fig. 0(a)] is 
separable. Moreover, regions in phase space that will not 
make contact with the SN-interface are in this geometry 
of measure zero, i.e. assumption A2 is valid. In sharp 
contrast to the discrepancy between Nqm and Abs, we 
observe that Aqm and Aebk agree almost perfectly in 
the whole energy interval < e < A (see Fig. CJ even 
for very elongated cavities where most orbits are non- 
retracing (Sy > Af). The criterion Sy w Af is found to 
be a good estimate where the BS quantization deviates 
from both the EBK quantization and the exact quantum 
spectrum. The case of the elongated A-billiard clearly 
demonstrates that the failure of the BS quantization in 
this case is due to the retracing assumption (Al), but not 
due to the failure of semiclassics (A3) itself. 
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C. Stretched non-separable billiard 

We focus now on billiards stretched in the direction 
of the SN-interface with D/lf » 1 and D/L » 1 
[Fig. Eft), c)]. For D/W > 1 the Andreev billiard becomes 
non-separable and an EBK quantization is no longer rig- 
orously justified. However, as we will show below, the 
predictions by the BS quantization can still be compared 
with an approximate EBK quantization. As the wave 
functions (Fig. (HJ) clearly indicate, the retracing approx- 
imation breaks down with large discrepancies between 
the particle (\u\ 2 ) and hole (\v\ 2 ) densities occurring. For 
this system, one obvious culprit is assumption A2. An ex- 
tended region of classical phase space (Fig.JUJ) remains de- 
coupled from the SN-interface. Consequently, the quan- 
tum DOS is associated, in part, with decoupled regions 
which are not represented at all by the BS approximation. 
It has previously been pointed out for other geometries 
(whispering gallery trajectories in circular billiards) that 
Andreev billiards indeed feature quantum states that are 
not associated with periodic orbits. 15 

A more detailed inspection reveals that a certain class 
of wave functions has the following features: (1) The 
probability density at the SN-interface is typically low 
- corresponding to decoupling from the superconductor. 
This suggests that the SN-interface effectively acts like 
a hard-wall boundary. (2) Integrating the probability 
density in the electron- and hole-part of the N-region 
typically shows a strong asymmetry. For the case de- 
picted in Fig.[S{c) we calculate a probability of 98.3% for 
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FIG. 6: (Color online) The DOS calculated quantum mechan- 
ically (red solid curve) and by the BS approximation (green 
dashed curve) for two different geometries with ratios of (a) 
L/D = 1 and (b) L/D = 6 (see insets) in units of the Weyl ap- 
proximation of the DOS per unit area, pw — rn e ft / (h 2 tt) . (c) 
error of the semiclassical prediction Sp [Eq. (|T5Jl] as a function 
of an energy- averaged r [Eq. JT1J] for 300 different ratios of 
L/D G [1, 20] . The solid blue line shows a smoothed average 
(30 adjacent points) of the recorded data. 
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FIG. 7: (Color online) Comparison of the exact quantum me- 
chanical state counting function A^qm (red staircase) with the 
semiclassical BS (green dashed line) and EBK (blue squares) 
approximations for a highly elongated cavity (L/D = 10), as 
shown in the lower inset. Quantum mechanical and EBK so- 
lutions are nearly indistinguishable. The upper inset shows a 
magnification. The energy for which 5y = Af is marked. 



the hole component and only 1.3% for the particle com- 
ponent. The coupling with the superconductor is very 
weak, resulting in a probability of 0.4% of finding the 
particle in the superconducting lead. These observations 
suggest that the particles in such non-retracing states are 
quasibound in either the electron- or the hole-space with 
only infrequent and short excursions to their respective 
mirror-space. As a consequence we may employ again an 
EBK quantization scheme to estimate the correspond- 
ing energy level of the Andreev state, however now in 
the space of the quasi-bound particle or hole. Exploit- 
ing the assumption of weak coupling to the SN-interface, 
we replace the SN-interface by a hard-wall (box) bound- 
ary. Applying this approach e.g. to the state depicted in 
Fig.JSJa), the EBK quantization 

_ 1 /n7r\ 2 1 /m7r\ 2 , N 

n ' m = 2 (t) + 2 (if) (16) 

implies for the density of the "isolated" electron state 
n = 6 maxima in x-direction and m = 205 maxima 
in ^/-direction. Inserting these quantum numbers yields 
an excitation energy of £box = 0.753 A, which is very 
close to the exact eigenenergy of the Andreev eigenstate 
of £qm = 0.760 A. Analogously, the state depicted in 
Fig. Etc) has quantum numbers n = 1, m = 41, corre- 
sponding to a hole excitation of £box = 0.038A, which 
compares well to the eigenenergy of the Andreev eigen- 
state of sqm = 0.040A. While the position of these An- 
dreev states can thus be explained fairly accurately by 
a box (or EBK) quantization, they are obviously not in- 
cluded in a standard BS approximation which considers 
only retracing electron-hole orbits. This illustrates the 
failure of assumption A2, while semiclassics still remains 
applicable. 
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FIG. 8: (Color online) Selected Andreev eigenstates of non- 
separable Andreev billiards stretched parallel to the SN- 
interface displaying a large discrepancy between electron and 
hole wavefunction patterns, (a) L/D = 10, e = 0.760A; 
(b) L/D = 10, e = 0.766A, k F = 21.bir/W; (c) L/D = 2, 
e = 0.040A, k F = 20.5tt/W; 
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FIG. 9: (Color online) (a) Poincare surface of section of a 
stretched non-separable Andreev billiard, taken at the right 
vertical side, opposing the SN-interface: W/D = 5, W = L. 
The light gray area marks regions of classical phase space 
coupled directly to the entrance lead, with one single reflec- 
tion at the right wall. Dark red areas show regions of classical 
phase space decoupled from the superconducting lead. Classi- 
cal trajectories corresponding to selected areas of phase space 
are shown in (b)-(d). 



(a) electron hole 




FIG. 10: (Color online) Selected Andreev eigenstates of sys- 
tems with (a) no diffractive corners (L/D = 1.6, W — D, 
k F = 20.5tt/W, e/A = 0.48); (b) one diffractive corner 
(D/L = 1.7, W = L, k F = 21.5tt/W, e/A = 0.72); (c) 
two diffractive corners (D/L = 1.7, W = L, k F = 20.5^/^, 
e/A = 0.82). 



For the present hard- wall TV-cavities, introducing non- 
separability implies simultaneously the introduction of 
diffractive edges and corners (see Fig. 0}. 36 In the pre- 
ceeding example, diffractive effects were present but for 
the special group of states considered of minor impor- 
tance. The latter was confirmed not only by the density 
distributions in the Andreev states (Fig. |SJ) but also by 
the accuracy of the EBK prediction. We can now en- 
hance diffractive effects by considering only moderately 
stretched cavities D/W > 1 with arbitrary position of 
the SN-interface. As shown in Fig. El we can thus real- 
ize Andreev billiards with zero corners (D/W — 1), one 
corner (D/W > 1, lower boundary of S and N aligned), 
and two corners (D/W > 1, arbitrary position of the 
SN-interface). Diffractive effects originate from potential 



variations on a length scale smaller than the de Broglie 
wavelength. In a semiclassical picture, diffractive scatter- 
ing gives rise to a non-deterministic, stochastic angular 
distribution. It certainly prevents the deterministic per- 
fect retracing of orbits. 

Diffractive scattering can be directly visualized in the 
wave function [Fig. irUT b.c)]. The density enhancement 
representing a scar of a classical trajectory hits the cor- 
ner of the SN-interface. As a result the trajectory "splits" 
into a bouncing-ball orbit and a broad angular distribu- 
tion. Significant differences between the particle (\u\ 2 ) 
and hole (\v\ 2 ) density in the retracing orbit are clearly 
visible. 

For a more quantitative assessment we compare again 
the quantum and the BS-results for the DOS [see 
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Fig. fTTT a.bl] and find, similarly to the case of the sep- 
arable geometry, that for larger D the mismatch Sp in 
the DOS increases. If we compare, however, how the er- 
ror Sp scales with the geometric retracing mismatch Sy 
[Fig. Illf c)]. we find that in the Andreev billiard with a 
diffractive corner the average mismatch is considerably 
larger than for the geometry without a sharp corner. 
The difference in Sp at fixed Sy can be thus attributed 
to diffraction. Note that diffractive scattering is an addi- 
tional mechanism for the failure of the retracing approx- 
imation, which is independent of the failure induced by 
long electron- hole orbits. 



o 0.5 e/A i o 0.5 s/A i 




FIG. 11: (Color online) The DOS calculated quantum me- 
chanically (red solid curve) and by the BS approximation 
(green dashed curve) for two different geometries with ratios 
of (a) L/D = 1 and (b) 1/8 (see insets), in units of the Weyl 
approximation of the DOS per unit area, pw = m e ft / (h 2 tt) . 
(c) Error of the retracing approximation Sp [Eq. (|15f> ] for 
transversely elongated cavities (red circles) as a function of 
r. The solid blue (dashed red) line shows the average error 
for longitudinally (transversely) elongated cavities (see also 
Fig. EJ. 

Additional evidence along these lines is presented in 
Fig. IT^l where the mismatch Sp between the BS approx- 
imation and the quantum results as a function of the 
geometrical aspect ratio of the normal-conducting region 
is compared for three geometries with either (a) zero, (b) 
one, or (c) two diffractive corners at the SN-interface. As 
expected, the mismatch increases (1) for larger aspect ra- 
tios and (2) at fixed aspect ratio for an increasing num- 
ber of diffractive corners. Further support comes from 
the observation that the mismatch Sp for the cavity with 
two diffractive corners is, to a large degree, independent 
of the position of the lead, as shown in Fig. 

Unlike other sources for the failure of the retracing ap- 
proximation discussed above, diffraction limits the valid- 
ity of standard semiclassical approximations itself, i.e. of 
assumption A3. Methods for describing diffractive effects 
semiclassically by including pseudopaths^LiS and non- 
deterministic scattering 3 ^ have been explored for scat- 



tering at open TV cavities. They play a crucial role in 
determining the S- matrix for such devices. The present 
results demonstrate that diffractive corrections come into 
play also for Andreev systems. 



IV. CONCLUSION 

We study Andreev billiards of rectangular shape and 
compare the quantum solutions of the Bogoliubov-de 
Gennes equation with a semiclassical Bohr-Sommerfeld 
approximation that assumes retracing electron-hole or- 
bits as a result of perfect Andreev reflection. Investi- 
gating the validity of this widely used semiclassical ap- 
proach, we identify three independent mechanisms which 
influence the density of bound Andreev states beyond the 
simple BS quantization: (i) long trajectories in the bil- 
liard typically magnify deviations between the electron 
and the hole part of orbits and cause a failure of the 
retracing approximation. This is of importance for the 
understanding of the failure of the BS approximation to 
reproduce the hard gap for chaotic N cavities. The lat- 
ter is controlled by the behavior of long paths for which, 
as we show, the retracing approximation breaks down, 
(ii) BS quantization fails when a subspace of the phase 
space is decoupled from the SN-interface and thus inac- 
cessible to retracing orbits, (iii) diffractive effects limit 
the validity of the semiclassical quantization. In the lat- 
ter case the failure is not just due to the assumption of 
retracing, but due to the failure of standard semiclassics 
itself. In this case, also alternative methods such as EBK 
quantization of separable or adiabatic quantization^ 3 ^ 
of (smooth) chaotic systems will fail. Corrections beyond 
standard semiclassics, e.g. inclusion of pseudo-paths 37,38 
of diffracted orbits will be needed. 




1 1.5 D/W 2 

FIG. 12: (Color online) Error of the retracing approximation 
Sp as a function of D/W for geometries with no (blue square, 
dashed line), one (red circle, solid line) and two (green cross, 
dashed line) diffractive corners. Different orientations of the 
cross refer to different positions of the superconducting lead. 
The lines show a linear root mean square fit of the data. 
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